Formation of regular spatial patterns in ratio-dependent predator-prey model driven 

by spatial colored-noise 
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Results are reported concerning the formation of spatial patterns in the two-species ratio- 
dependent predator-prey model driven by spatial colored-noise. The results show that there is 
a critical value with respect to the intensity of spatial noise for this system when the parameters 
are in the Turing space, above which the regular spatial patterns appear in two dimensions, but 
under which there are not regular spatial patterns produced. In particular, we investigate in two- 
dimensional space the formation of regular spatial patterns with the spatial noise added in the side 
and the center of the simulation domain, respectively. 

PACS numbers: 87.23.Cc, 05.40.-a, 87.10. +c 



I. INTRODUCTION 

The influence of noise on nonlinear systems is the sub- 
ject of intense experimental and theoretical investiga- 
tions. The most well-known phenomenon is noise induces 
transition [1] and stochastic resonance [2, 3], both show- 
ing the possibility to transform noise into order. 

In this paper, the problem of the predator-prey dynam- 
ical system [4, 5, 6, 7, 8] is revisited and the results are 
presented concerning the noise aspects of the formation 
of spatial patterns in two dimensions. To our knowl- 
edge, the predator-prey dynamics in the absence of noise 
have been studied extensively using mathematical mod- 
els. A common feature of these models is the prediction 
that the populations can cycle for some parameter sets: 
i.e. the populations of predators and prey do not set- 
tic to constant values, but rather, oscillate periodically 
in time. There are considerable data from field studies 
and laboratory experiments to support the existence of 
such population cycles. If a small group of predators are 
introduced into an spatial uniform population of prey, 
the predators will tend to invade the prey, leaving be- 
hind a mixture of predators and prey. In Jonathan A. 
Sherratt, Mark Lewis, Barry Eagan and Andrew Fowler's 
works [9, 10], they have studied the behavior of such in- 
vasions for cyclic populations. The results are somewhat 
surprising: the invasion leaves behind spatiotcmporal os- 
cillations falling into either periodic traveling waves or 
spatiotemporal irregularity; mixed behavior also occurs. 

One of the key issues in ecology is how environmental 
fluctuation and species interaction determine variability 
of the population density [11, 12, 13]. The other topic 
is that dynamic patterns in two-dimensional space have 
recently been introduced into ecology [14, 15, 16, 17]. 
Ecologists have mainly been interested in the dynami- 
cal consequences of population interactions, often ignor- 
ing environmental variability altogether in space. How- 
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ever, the essential role of environmental fluctuations has 
recently been recognized in theoretical ecology. Noise- 
induced effects on population dynamics have been sub- 
ject to intense theoretical investigations [8, 18, 19, 20, 21, 
22, 23, 24, 25, 26]. Moreover, ecological investigations 
suggest that population dynamics is sensitive to noise. 
In spite of the obvious significance of this circumstance, 
the role of non-equilibrium spatial fluctuations (spatial 
noise) of environmental has not been investigated much 
in the context of ecosystems. 

Recently, the general N-species Lotka- Voltcrra and 
ratio-dependent predator-prey dynamical systems are de- 
scribed as stochastic models by Romi Mankin, Ako Sauga 
and et al. i.e. add the parameter noise term in the classi- 
cal models, where the carrying capacity of the population 
is taken into account as dichotomous noise. For former 
model, the result of the study has shown that the noise 
induces discontinuous transitions from a stable state to 
an instable one in Refs. [21, 22]. The latter model [8] has 
shown that the phenomena of colored- noise-induced Hopf 
bifurcation and corresponding reentrant transitions can 
appear in ecosystems. But both models are studied in 
nonspatial. From the Refs. [9, 27, 28] we know that the 
space plays a crucial role in the ecosystem. In Rcf. [28], 
the results have shown that the Turing patterns occur in 
the ratio-dependent prcdator-prcy models. Thus, natu- 
rally arise a problem that how the spatial noise influences 
the formation of spatial patterns or species invasion in the 
ratio-dependent predator prey dynamical systems? Here 
we address how the spatial noise affects the formation of 
the spatial patterns on the coexistence equilibrium in the 
reaction-diffusion predator-prey dynamical systems. 



II. MODEL 

A. The ratio-dependent model with diffusion 

We consider the reaction-diffusion model for predator- 
prey interactions. At any point (x, y) at time T, the 
dynamics of prey (N(x,y;T)) and predator (P(x,y;T)) 
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populations are given by a reaction-diffusion model with 
logistic growth of the prey and a type II (or Michaelis- 
Menten-type) function response of the predator 
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where V 2 is the Laplacian operator in Cartesian coordi- 
nates. The parameters r, K, fi and e denote the intrinsic 
growth rate, carrying capacity of the prey, the death rate 
of the predator and the yield coefficient of prey to preda- 
tor, respectively. The constant D\ and D 2 are prey and 
predator diffusion coefficients, respectively. 

The model (1) can be simplified by introducing the 
dimensionless variables, n = N/K, p = Pk/K, t = Tr, 
x = x/y/Di/r, and y = y/yjD\/r. The units of the 
dimensional variables will be scaled. Thus, in the terms 
of these dimensionless variables the model is simplified 
to 

| = (l-n)n-^ + VV (2a) 
at n + p 

^=eP^- VP + dW 2 p, (2b) 
at n + p 

where four new dimensionless parameters naturally arise: 
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(3) 



From the Ref. [4], we know that the ratio-dependent 
predator-prey dynamical system (2) has a homogeneous 
coexistence point, (n s ,p s ), i.e. the coexistence point is 
given by 



n s = 1 - {r)/s)(A),p s = (A - l)n s 



(4) 



where A = e/3/rj. This coexistence point is biologically 
meaningful (n s > 0,p s > 0) only if r\ > e(fi — 1) and 
rj < ef3 

Since we focus on the spatial noise affects the formation 
of the spatial patterns in the ratio-dependent predator- 
prey system, we don't show the mathematical property 
here. 



B. Spatial noise in predator-prey model 

Environmental fluctuations and other factors (see the 
discussion) are important components in an ecosystem. 
In Refs. [29, 30], May have pointed out the fact that 
because of environmental fluctuations, the carrying ca- 
pacities, death rates, birth rates, and other parameters 
involved with the systems exhibit random fluctuations to 
a greater or lesser extent. Because a multinoise fluctua- 
tion stochastic system is difficult to consider analytically, 



it is generally reasonable to confine only one (or two) pa- 
rameter of the system to fluctuate. In our paper, a ran- 
dom interaction with the environment (such as climate, 
disease, etc) in some domain is taken into account by in- 
troducing a random spatial noise (£(x,y;t)) to the prey. 
Henceforth, the dynamic system (2) is used to numerical 
simulation. In comparison with Refs. [8, 21, 22, 23], here 
the difference is that we take into account the spatial 
noise in two dimensions. For the spatial noise we choose 
the colored noise £(x, y;t), which accounts for density 
fluctuations at point (x,y). The Ornstein-Uhlenbeck pro- 
cess was used to calculate the colored-noise (Appendix). 
To characterize the spatial noise, we introduce parame- 
ter, D, the intensity of the noise. The £(x, y;t) has the 
following properties: 



<£(x,y;t)> = 



and 



{fa, y; t)t(x>, y'; t')) = D 2 S(x - x')S(y - y^-^. 

Obviously, r is the noise correlation time. Although 
the dynamical system (2) absent of random spatial noise 
term is well understood, there is little analytic theory 
for the system in two-dimensional space. The equations 
are integrated numerically with an implicit scheme us- 
ing 200 x 200 an 50 x 200 spatial grid sites, respectively. 
For nonlinear equations, the discretization introduced by 
numerical methods may generate spurious results. To 
test this possibility, the resolution of the simulation re- 
sult was increased with space and time, and the system 
is integrated with a different scheme, a finite difference 
method. In all case the same qualitative results have been 
obtained. We assume zero flux boundary conditions for 
all t in the simulation. 



III. RESULTS 

We study the effect of the spatial noise to the forma- 
tion of the spatial patterns. The results presented here 
are based on numerical analysis for a set of parameters, 
= 2, £ = 0.5, 7] = 0.6, d = 10 [31], Ah = 0.5 (space 
scales), and At = 0.005 (time step), chosen to obtain 
large-scale spatial patterns in the absence of spatial noise 
(see Fig. 1). Moreover, from Turing Structures theory we 
know that these parameters must be within the Turing 
space. In fact, the authors [27] have validated that these 
parameters are in the Turing space. 

Figure 1 illustrates the irregular spatial patterns be- 
havior of prey numbers in two dimensions absent of the 
spatial noise. From this figure one can see that the spa- 
tial patterns are irregular in two dimensions when they 
reach the steady state. 

The main results in our paper show the effect of spatial 
noise with respect to formation of the spatial patterns. In 
order to ensure the facticity of the result, we have studied 
the system (2) in two cases driven by the spatial noise. 
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FIG. 1: (Color online) Spatial patterns in two di- 
mensions absent of noise (in this case, the qualita- 
tive results are the same as the Ref. [27]). The nu- 
merical integration of the dynamical system described 
by Eqs. (2) has been performed in a 200 x 200 lat- 
tices. (A) t = and (B) t = 550. The right in- 
set shows the color corresponding to the value of prey. 




FIG. 2: (Color online) Spatial patterns in two dimen- 
sions driven by spatial noise. The numerical integra- 
tion of the dynamical system described by Eqs. (2) has 
been performed in 50 x 200 and 200 x 200 lattices, 
respectively. (A) The spatial noise is added at the 
left-hand side of the 50 x 200 lattice; (B) The spatial 
noise is added at the center of the 200 x 200 lattices. 



First, we introduce the spatial noise (here the spatial 
noise is perturbation on the equilibrium value) at the 
left of the 50 x 100 spatial grid sites. The domain of the 
spatial noise is set by 50 x 3. The ecological significance 
of the spatial noise can be explained for the invasion of a 
prey population by predators in the space (see Fig. 2(A)). 
Second, we introduce the spatial noise (here the spatial 
noise is the same as that in first) at the center of the 
200 x 200 spatial grid sites (see Fig. 2(B)). 

Figure 2 illustrates the formation of the regular spa- 
tial patterns behavior of the prey numbers in two dimen- 
sions after the spatial patterns reach the steady state, 
where the initially condition (at t = 0) add random and 
nonuniform small perturbations to the equilibrium val- 
ues. The maximum amplitude of the perturbations is 
less than 0.0005 and the noise correlation time r = 1. 
The intensity of the spatial noise is equal to 0.05. Com- 
paring Fig. 1 with Fig. 2, one can easily find that some 
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FIG. 3: This graph illustrates the regular pat- 
terns driven by the spatial noise from other views. 

characteristics are the same (here mainly relate to the 
wavelength). But Fig. 2(A) and Fig. 2(B) both appear 
periodical structure patterns in two dimensions (called 
the spatial-period-2-structure) . Additionally, these spa- 
tial patterns are regular in two dimensions. The only dif- 
ference is that In Fig. 2(A) the spatial patterns are stripe 
but in Fig. 2(B) target when the system reach steady 
state. 

If we change the direction of observing angle, the space 
vs density plot can be seen more clearly. The space plot 
of the variable n is shown in Fig. 3, where the horizontal 
axis represents the horizontal spatial line crossing the 
center point. Fig. 3 corresponds to Fig. 2(B). 

Finally, we increase the intensity of the spatial noise 
in the system, and the results show that for every ini- 
tial condition there is a critical value of the intensity of 
noise, D c . Here the critical value D c ss 0.02 with fore- 
going parameter. Above the critical value the regular 
spatial patterns appear and maintain in two dimensions. 
However, under the critical value, there doesn't exist reg- 
ular spatial patterns even though the simulation runs a 
long time. The qualitative results are the same under 
different D values (above the critical value), so we only 
show typical results as in Fig. 2 when D = 0.05. Under 
the critical value, the result is not shown in our paper 
since the spatial patterns are similar to Fig. 1. 

IV. CONCLUSION AND DISCUSSION 

In our paper we have considered the effect of the spa- 
tial noise with respect to the spatial patterns in the 
ratio-dependent predator-prey system in two dimensions. 
The results show that there exits a critical value of the 
intensity of the spatial noise (D c ) for the system (2) 
when the parameters are in the Turing space. Above 
the critical value the regular spatial patterns appear in 
two-dimensional space, but under the critical value there 
doesn't exit regular spatial patterns. Our key objective 
is that the regular spatial patterns can appear in ratio- 
dependent predator-prey model driven by spatial noise. 
Here, we use the intensity of the noise to characterize the 
spatial noise, one can also use other measures to charac- 
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terize the spatial noise (e.g signal to noise ratio). 

Our results differ in several ways from the previous 
studies of predator-prey equations in ecology [8, 20, 21, 
22, 23, 32]. Previous studies have shown that the diffu- 
sion can induce instability and the parameter noise can 
induce Hopf-bifurcation, discontinuous transition and so 
on. But here we investigate the spatial noise in the sys- 
tem (2) in two dimensions. Our results suggest a direc- 
tion for further research, which is how to detect the spa- 
tial fluctuation of the complex dynamics in natural envi- 
ronments. In addition, two species predator- prey systems 
exhibit Turing spatial patterns only under special condi- 
tions. In Ref. [27], the study suggests that any strat- 
egy that involves interference between predators should 
enhance patterns formation. Furthermore, the process 
rates (prey birth rate, predator death rate, predator at- 
tack rate) and the mobility of population size can also 
be determining factors in the formation of biological pat- 
terns. 



Appendix 

The colored noise satisfies 



It can be shown that the colored noise £(f) my be calcu- 
lated form the following equation 



dt 



1 L>V2t 



(6) 



where 77(f) is a Gaussian white noise source with 

(v(t)) = 0, 
(r?(f)7?(f)) = 5(t-t>). 

Equation (6) describes an Ornstein-Uhlenbeck pro- 
cess [33], and it can be shown that the solution to this 
equation satisfies Eqs. (5). 
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